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We work in the framework of an inflationary cold dark matter universe with cosmological con- 
stant, in which the cosmological inhomogeneities are considered as gravitational lenses for the CMB 
photons. This lensing deforms the angular distribution of the CMB maps in such a way that the 
induced deformations are not Gaussian. Our main goal is the estimation of the deviations with 
respect to Gaussianity appeared in the distribution of deformations. In the new approach used in 
this paper, matter is evolved with a particle-mesh N-body code and, then, an useful ray-tracing 
technique designed to calculate the correlations of the lens deformations induced by nonlinear struc- 
tures is applied. Our approach is described in detail and tested. Various correlations are estimated 
at an appropriate angular scale. The resulting values point out both deviations with respect to 
Gaussian statistics and a low level of correlation in the lens deformations. 

PACS numbers: 98.70.Vc, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

In the absence of a reionization, the photons of the 
Cosmic Microwave Background (CMB) are not scattered 
from decoupling to present time, except inside galaxy 
clusters (Sunyaev-Zel'dovich effect); nevertheless, the 
propagation direction of these photons changes due to the 
gravitational action of the cosmological inhomogeneities 
(lensing) . 

The so-called primary CMB anisotropics were pro- 
duced by linear structures at high redshifts. The cor- 
responding temperature contrast, A^, = (AT/r)p, is a 
statistical field which can be expanded in spherical har- 
monics as follows: = 'J2ti=-e '^^rnYerm where the aem 
quantities are statistically independent variables depend- 
ing on the sky realization we are expanding. If many sky 
realizations are averaged, the resulting a^„i coefficients 
have zero means and variances (| aim |^} = Ci. Since Aj, 
is a homogeneous and isotropic statistical field, quanti- 
ties (I aim 1^) do not depend on m. In the model under 
consideration, scalar energy density fluctuations are ini- 
tially Gaussian, and they remain Gaussian during linear 
evolution, by this reason, the primary CMB anisotropics 
-which were produced by linear inhomogeneities- also are 
Gaussian. 

Gaussian primary anisotropy appears to be superim- 
posed to secondary anisotropics generated well after de- 
coupling, and also to the anisotropy of the microwave 
radiation emitted by both our galaxy and some ex- 
tragalactic sources (contaminant foregrounds). These 
foregrounds are sub-dominant non Gaussian components 
contributing to the total observable signal at microwave 
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wavelengths. In the model we are considering here, there 
are various secondary anisotropics. A few comments 
about them are worthwhile: the late integrated Sachs- 
Wolfe anisotropy was produced by quasi-linear inhomo- 
geneities Ij and, consequently, deviations with respect to 
Gaussianity are not expected to be very important; an- 
other secondary gravitational anisotropy was produced 
by strongly nonlinear structures as galaxy clusters and 
substructures (Rees-Sciama effect), this component is 
sub-dominant and non Gaussian, and its deviations with 
respect to Gaussianity were studied by various authors 
0, S ; finally, the Sunyaev-Zeldovich anisotropy was 
produced by hot plasma in galaxy clusters, and it is also 
non Gaussian 

The main goal of this paper is the analysis of the 
non Gaussianity caused by lensing, and such an analysis 
can be performed using the operative ray-tracing proce- 
dure designed in next sections. In the absence of lens- 
ing, there are directional frequency shifts which cause an 
anisotropic temperature distribution in the CMB. The 
lens effect does no cause frequency shifts, but it pro- 
duces angular deviations of the propagation directions. 
This means that the lens effect does not produce any al- 
teration of the CMB temperature. It changes the prop- 
agation direction of the CMB photons and, then, the 
angular distribution of the CMB temperature changes 
accordingly. In other words, the lens effect deforms the 
maps of the CMB temperature which appear as a result 
of pure frequency shifts. Any numerical estimate of the 
alterations of the CMB maps produced by lensing (in- 
cluding non Gausianity generation) involves three steps: 
(i) maps of the CMB temperature distribution are built 
up for the model under consideration (without any lens- 
ing), (ii) these maps are deformed taking into account 
the deviations of the propagation directions produced by 
lensing, and (iii) the resulting deformations of unlensed 



2 



maps (hereafter called either lens deformations or lens 
effect) are analyzed. Let us now consider each of these 
steps in more detail. 

In order to built up unlensed maps of the CMB, we 
assume that, in the absence of lensing, the Gaussian pri- 
mary anisotropy dominates and, consequently, the maps 
to be deformed can be built up using the angular power 
spectrum of this primary anisotropy, which has been cal- 
culated using CMBFAST 6] for the ACDM model under 
consideration. As it is discussed below, only small Gaus- 
sian maps having sizes of a few degrees and resolutions 
of the order of one arc-minute are necessary and they 
are built up using the Fourier Transform (FT) and the 
mentioned spectrum {Ci quantities). See Q and refer- 
ences cited therein for details. The resulting maps have 
spots with different sizes and amplitudes (angular struc- 
ture). Spots with angular size 7t/£ have a mean am- 

1/2 

plitude proportional to ; hence, spots with angu- 
lar size close to one degree {£ ~ 200) have large am- 
plitudes, whereas spots corresponding to one arcminute 
{£ ~ 10000) have negligible amplitudes as a result of the 
smallness of Cioooo- For the spectrum under considera- 
tion, significant spots have sizes of various arcminutes. 

Cosmological inhomogeneities produce lensing and 
this effect is associated to the deviation field S. The 
anisotropy observed in the n direction is the primary 
anisotropy corresponding to the direction uq in the ab- 
sence of lensing, where ft = fto — S; hence, 

A{n) = A,{no) = A^{n + S). (1) 

The unit vectors n and no point towards two points of 
the last scattering surface, and the deviation field 6 gives 
the angular excursion on this surface due to lensing. We 
must estimate the deviation field i5 (see below) to get 
the deformed observable map A. If the 6 deviations are 
smaller than the size of the smallest significant spots (var- 
ious minutes), Eq. can be expanded in S powers to 
get 

Ain)^A,{n) + ^-S (2) 

We see that map deformations, namely, the differences 
in the temperature contrasts corresponding to associate 
directions (second term of the r.h.s.) are the product of 
two statistical fields: A^, and 6 and, even if these field 
are Gaussian and statistically independent, the resulting 
product would be non Gaussian. 

Several authors have studied angular correlations in 
lensed CMB maps. In the most theoretical papers on 
this subject, some correlations have been estimated with- 
out using simulations; the approach used in these papers 
only applies for those correlations which can be writ- 
ten in terms of the time varying form of the matter 
power spectrum. In the most basic of these papers Q, 
the mentioned spectrum was modelled beyond the linear 
regime and the second order correlations of lensed CMB 



maps were estimated. A similar method was used by 
Bernardeau to estimate fourth order correlations cor- 
responding to various sets of four directions (with distinct 
configurations); furthermore, this last author pointed out 
that third order correlations vanish. Other methods are 
necessary to estimate correlations of order n for n > 4, 
and also to simulate CMB maps with the true statistics 
induced by lensing, namely, maps having the true correla- 
tions at any order (not only at second and fourth orders) . 
The most powerful technique to study lens distortions 
in CMB maps is the use of ray-tracing through N-body 
simulations. Two methods to app ly this tecnique are de- 
scribed in references |l3 and |0| (other relevant papers 
in this field are cited in these two references). Another 
new method is proposed here. We agree with White & 

Hu who stated that: On subdegree scales, a full 

description of week lensing therefore requires numerical 
simulations, the most natural being N-body simulations 
; in fact, using ray-tracing through N-body simu- 
lations, CMB photons are deviated by fully nonlinear 
structures (as galaxy clusters) which move in the simula- 
tion box. Nowadays, only N-body simulations lead to a 
proper description of a distribution of strongly nonlinear 
structures. Unfortunately, these simulations constraint 
us to work in a periodic universe, and relevant problems 
associated to periodicity must be solved. Each method - 
designed to use ray-tracing through N-body simulations- 
corresponds to a different way for preventing periodicity 
effects. In this paper, one of these methods is described 
and tested. 

Hereafter, quantity h is the reduced Hubble constant 
h — 10~^Ho, where Hq is the Hubble constant in units 
of Kms'^^Mpc~^ , the density parameters corresponding 
to baryons, dark matter, and vacuum, are fib, i^d and 
r^A, respectively, the total density parameter is f2b + 
fld + ^A, and the matter density parameter is Qm = + 
fid- In the flat inflationary universe under consideration, 
the above parameters take on the following values: h — 
0.65, Qb = 0.05, nd = 0.25, and flA = 0.7 and, then, 
according to Eke et al. 12], the power spectrum of scalar 
energy density perturbations must be normalized with 
as — 0.93 in order to have cluster abundances compatible 
with observations. Units are chosen in such a way that 
c = SttG — 1, where c is the speed of light and G the 
gravitation constant. Whatever quantity "A" may be, 
Aj^ and stand for the A values on the last scattering 
surface and at present time, respectively. The scale factor 
is a{t), where t is the cosmological time, and its present 
value, ao, is assumed to be unity, which is always possible 
in fiat universes. 



II. FORMALISM 

The CMB photons move on null geodesies and the line 
element is: 

ds2 ^ _^ 2(t))dt^ + (1 - 2(j))a^Sijdx'dx^ (3) 
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where function (f) satisfies the equation: 



1 



(4) 



and pmo = Peril is the background energy density for 
matter. On account of this last equation, function can 
be interpreted as the pecuhar gravitational potential cre- 
ated by the cosmological structures. Equations Q and 
Q are valid for linear inhomogeneities located well inside 
the horizon and also for nonlinear structures (potential 
approximation), see [l^. 

The deviation field 5 is given by the following integral 



6 = -2 



(5) 



where V^c/) = — r? A n A V0 is the transverse gradient of 
the potential, and W{X) — {X^ — X)/X^. The variable A 
is 



X{a) = 



db 



(a„0&+r!A54)l/2 



(6) 



The integral in the r.h.s. of Eq. ||SJ) is to be evaluated 
along the background null geodesies. In our flat back- 
ground, the equations of the null geodesies passing by 
point are: 



X{a)n' 



(7) 



In next section, spherical clusters with Navarro-Frenk- 
White (NEW) density profiles are considered for 
qualitative analysis. These profiles appear in N-body 
simulations of dark matter halos and they have the form: 



p(r) 



Sc PmO 



(8) 



where pmO = Pcrit^mo, and pcrit = iH^ is the critical 
density of the universe. The parameters 5c and Vs are 
related to both the concentration parameter c„p„- and 
the mass M2oo- This mass is that contained inside a 
sphere whose mean density is 200. The radius of this 
sphere is denoted r2oo- The following relations hold: 



Ts = (3M2oo/8007rp„,o)'/Vc 



(9) 



CjvFw = 6-8, cluster C 



[T^. Fortunately, our qualitative analysis -presented in 
next section- does not require any detailed z-dependent 
profile. 

In Sec. Ill, three clusters -which are hereafter called 
Cj, Cjj, and C^j^- are considered. Cluster is a rich 
one having M200 = 1000 Mg (where Mg = 10^^ M©) and 
is a standard cluster having a few 

= 7.7, 

and cluster C^^^ is a galaxy group having a few tens of 
galaxies with M200 = 25 Mg and c^^^^^, = 12.5. Since the 
NFW profile diverges as physical radius r tends to zero, 
an uniform core of radius Tc has been assumed and, con- 
sequently, the density profile has the form ||SJ) for r > Tc 
and a constant value for r < re- Both profiles match con- 
tinuously at r = For this mass distribution, function 
dip/ dr has the following form: 



cents of galaxies with M200 = 250 Mg and c 



1 



p{rc)r 



dr 6' 

for r > Vc, where p{rc) is given by Eq. (jSJ, and 



(11) 



(12) 



for r < rc, where A = (l/6)rj?p(rc), D ^ 6c Pmo/2, and 
E{r) = ln[l + (r/rs)] + [I + {r/rs)]~^. Core radius of 
0.2 Mpc, 0.1 Mpc, and 0.05 Mpc are assigned to clus- 
ters Cj , Cjj , and Cj^j , respectively. In next section, our 
conclusions are probed to be highly independent on this 
assignation. 

Spherical linear inhomogeneities evolving inside the ef- 
fective horizon are also considered in Sec. Ill for qual- 
itative discussion. These spherical regions are assumed 
to be uniform and they are characterized by their co- 
moving radius i?^ . This radius fixes the associated 
comoving scale k^. = 7r/2i?p and, then, the present 
density contrast is fixed using the power spectrum of 
the energy density perturbations -in the model under 
consideration- which gives the average amplitude of the 
Fourier mode . The density contrast evolves propor- 
tional to the growing mode D(t), which is given in .IT^ . 
Three of these linear inhomogeneities will be considered, 
they correspond to = 250 Mpc, = 50 Mpc, and 
R^, = 12.5 Mpc ^ 8h^^ Mpc, and they are denoted Lj, 
, and L^^j , respectively. The derivative of the gravi- 
tational potential with respect to the physical radius r is 



200ct,„ , , dr 6 

(10) 



^ ^lip- Pmo)r (13) 



ln[l 4- Cr^piv ) /(I + (^NFw ) 



for r < aR^ , and 



c 



The profile |(SJ) remains almost unchanged from virializa- 

tion (z ~ 1) to present time. Some authors have tried to d(f> a'^R^ 

assume the NFW profile with variable values of the in- d-f ~ j-2 ^j. 

volved parameters to go beyond virialization -until z ~ 5 

in reference |l5j | -. Other authors use generalized profiles for r > aR^ . 



(aRa) (14) 
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III. ANGULAR SCALES, BOXES AND 
RESOLUTION 

Recently, various groups have used simulations of 
structure formation to estimate the effects of lensing on 
the CMB. An important problem is that there are many 
scales producing lensing and that, according to White & 
Hu |Ij , simulating the full range of scales implied is cur- 
rently a practical impossibility. These authors proposed 
the tiling method to circumvent this problem and also the 
problems with periodicity (see Sec. I); this method em- 
ploys many independent simulations with different boxes 
and resolutions to tile the photon trajectories. White & 
Hu 1^ explained the advantages of their method with 
respect to the most traditional one based on plane pro- 
jections, which has been extensively used in the literature 
(see ^3 and references cited therein). The method pro- 
posed here does not require many independent simula- 
tions -but only one- and, consequently, it is simpler and 
numerically faster than the tiling one. Furthermore, it 
does not use plane projections whose potential problems 
were pointed out in pJ| . In this section, we are con- 
cerned with the scales relevant for CMB lensing, whereas 
periodicty is considered in next sections. 

We are interested in the deviations with respect to 
Gaussianity produced by lensing These deviations 
can be measured by angular correlations corresponding 
to three and more directions. Given an angular scale, the 
correlations are produced by a set of linear and nonlinear 
cosmological structures and, a detailed study of the type 
of structures contributing to these correlations is neces- 
sary in order to design our N-body simulations, which 
should include -in the same box- all the structures re- 
sponsible for the effects we are looking for. In order to 
estimate scales, we can consider the evolving linear in- 
homogeneities Lj, L^^, and i^jj, and the nonevolving 
NFW density profiles Cj , C^^ , and C^^j . Cluster evolu- 
tion would not affect our qualitative estimates for small 
redshifts z < 5, although it could be crucial in other con- 
texts. In short, our study is based on six appropriate 
structures, three linear ones and three galaxy clusters. 

First of all, Eq. and the potentials given in Sec. II 
have been used to calculate the lens deviations produced 
by the six selected structures. Each linear structure is 
placed at redshifts 2, 10 and 100, whereas each cluster is 
located at redshifts 0.5, 2 and 5 (cluster Cn is also placed 
at z = 0.2 and z = 0.1). The deviation 6 is calculated 
for each observation direction, which is characterized by 
the angle 6 formed by the line of sight and the direction 
pointing towards the structure center. Results are pre- 
sented in Fig. n where we see that: (i) the deviations 
grow as the redshift decreases, and the growing is small 
for z ^ 0.5 (compare the curves of the left middle panel 
of Fig. n corresponding to redshifts 0.5, 0.2 and 0.1), 
(ii) the deviations increase as the inhomogeneity size in- 
creases (this size is fixed by M200 for clusters and by i?^ 
for linear structures), (iii) small linear structures {L^^j) 
produces very small deviations whatever their location 



may be, (iv) small clusters (or groups, C^^j) produce a 
maximum deviation which is about 10 % of the maxi- 
mum deviation produced by cluster Cj (located at the 
same redshift), and (v) the most important deviations 
are produced by massive clusters and by very extended 
linear inhomogeneities. 

We are interested in the angular correlations of some 
CMB maps. Given a map of the variable C and m direc- 
tions, these correlations are defined as follows: 

C™ = (C("i)C("2)---C(^™)) , (15) 

where the average is over many realizations of the CMB 
sky. Some of these averages are estimated below for 
two, three and four directions in maps of both primary 
anisotropy and lens deformations. For m = 4, the cho- 
sen directions draw the vertices of a square on the last 
scattering surface with fii ■ n2 = cos a, and for m — S, 
they point towards the vertices of an isosceles rectan- 
gle triangle with ni • n2 = cos a and ni ■ 713 — cos a. 
The relevance of the deviations presented in the panels 
of Fig. ^depends on the angular scale, a, chosen for cor- 
relation computations. Suppose the scale a it = -k/o) 
and m — 2. From the primary temperatures correspond- 
ing to many pairs of directions (ni, 712) forming angle 
a, the average (A^ (?ii)Ap (712)) can be calculated and 
the result measures correlations in the absence of lens- 
ing. According to Eq. (Q, after lensing, the correla- 
tions are (A(rri)A(?i^)) = (A^ (n^i) A^ (7102)} and, con- 
sequently, the correlations would change -due to lensing- 
if, as a result of deviations, the directions riQi and no2 
form an angle ao different enough from a; therefore, in 
order to see if a certain spherical inhomogeneity located 
at redshift z can contribute to the correlation at the an- 
gular scale a, we may consider two directions pointing 
towards R and S, two points of a certain diameter of the 
spherical structure; the first (second) of these directions 
forms angle 6 [6 + a ) with the direction of the sym- 
metry center, and the angle O = 9 + a/2 corresponds 
to the middle point of the segment RS\ then, the dif- 
ferences d{Q) = 6{0 + a) — 6{9) measure the deforma- 
tions of the angle formed by the two chosen directions. 
The top panels of Fig. |21 give these differences for the 
structures C^, C^^, and Cj^^ (at z=0.5), and the bottom 
panels for Lj and Lj, (at z=2). The correlation scale 
is a = 13.5' = 800) in the left panels and a = 6.75' 
~ 1600) in the right ones. From the analysis of this 
Figure it follows that: (a) big close clusters of type Cj 
produce the largest deformations reaching the maximum 
value ~ 1.68' (- 1.98') for a = 13.5' (a = 6.75'), (b) for 
a = 13.5', inhomogeneities like C^^-, Cj^^, Lj, and L^^ 
produce maximum deformations which are about 37.5%, 
7.7%, 33% and 6.8% of the maximum ones corresponding 
to clusters (~ 1.68') and, (c) for a — 6.75', the max- 
imum deformations produced by Cj, , Cjjj , Lj , and L^^ 
structures are about 36.9%, 7.1%, 13.7% and 2.9% of the 
maximum deformation appearing in case 1.98'). 
Points (b) and (c) lead to the conclusion that the de- 
formations produced by big linear structures are more 
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FIG. 1: Lens deviation (5 as a function of the angle 6 between the line of sight and the direction pointing towards the center of 
a spherical inhomogeneity. Top, middle, and bottom right (left) panels correspond to the linear structures L^, L^j, and Ljjj 
(clusters Cj, Cjj, and Cjjj) defined in the text. In right (left) panel, solid, dotted, and dashed lines correspond to redshifts 2, 
10, and 100 (0.5, 2, and 5), respectively. In the middle left panel, dotted-dashed (triple-dotted-dashed) line shows results for 
redshift 0.2 (0.1). 
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FIG. 2: Lens deformation, d, of the angle a formed by two directions pointing towards the same diameter of a spherical 
structure vs. the angle Q defined in the text. This angle fixes the relative position of these directions with respect to the line of 
sight of the symmetry center. Solid lines of the top (bottom) panels correspond to cluster at z = 0.5 (the inhomogeneity 
at z = 2), whereas dotted lines correspond to the same redshifts but different structures Cj^ and Lj^ in all the panels. Dashed 
hnes of the top panels are associated to cluster C^^^ at z = 0.5. Left (right) panels give deviations for a = 13.5' (a = 6.75'). 



important for a — 13.5', but they are not negligible in 
the case a = 6.75'. The deformations corresponding to 
structures like L^^^ are too small and they have not been 
presented in Fig. |2 

From the above discussion, it follows that -for the 
scales under consideration- the most relevant structures 
are clusters of type Cj and Cjj at low redshifts, and 
big linear inhomogeneities with > 50 Mpc located 
at small redshifts. This conclusion is important. It tell 
us that simulations which do not resolve galaxy groups 
and small clusters can lead to good estimations, and also 
that the box size can be chosen to prevent the existence 
of extended linear inhomogeneities contributing to lens- 
ing. Boxes with sizes of 128 Mpc should be appropri- 
ated, sizes of 256 Mpc could be acceptable (although 



a part of the lensing could be due to linear structures 
with i?p > 50 Mpc, see Sec. V) and, finally, boxes of 
512 Mpc and greater could contain an extended part of 
a structure giving important contributions to lensing. 
The idea is that the lens effect produced by big linear 
inhomogeneities can be studied without simulations, see 
references in whereas the effect of any other signif- 
icant structure (C^ and C^^ clusters) can be included in 
simulations for appropriate boxes (between 128 Mpc and 
256 Mpc size) . Let us now analyze if results from Figs, 
n and 121 arc robust against variations of the core radius 

Imagine one of the above three clusters located at red- 
shift z. It produces a significant deviation of the CMB 
photons if and only if these photons cross the region 
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FIG. 3; Function 2W{\)V±(f) x lO" -involved in Eq. Q- 
v.s. the comoving distance A — Ac to the cluster center for 
the observation angle 6 = 2.5'. All clusters have been located 
at 2 = 0.5. Each pair of associated solid and dotted lines 
corresponds to the cluster indicated inside the panel. Solid 
(dotted) line has been obtained from the greatest (smallest) 
values assigned to the core radius in the text. 



where the function to be integrated in Eq. |(SJ (propor- 
tional to V±(j>) is not neghgible; hence, if the dependence 
of this function on the core radius is weak, lens deviations 
(integral in Eq. |(SJ) also depend weakly on r^, and pre- 
vious results based on Figs, ^andl^jare robust. In order 
to study the mentioned function for clusters Cj , C^j , and 
Cjjj, these structures are located at redshift 0.5 (where 
cluster evolution is not expected to be important) and, 
then, the function d6{9)/dX = 2W{X)^±(j) involved in 
Eq. is calculated at points of the null geodesic cor- 
responding to the angle 9 = 2.5'. Results are presented 
in Fig. 13 where function d6{2.5)/dX is given (with arbi- 
trary normalization) in terms of the comoving distance 
to the cluster center A — Ac, where A is the comoving dis- 
tance from the observer to an arbitrary point of the null 
geodesic, and Ac is the A value corresponding to the clus- 
ter center. Continuous lines display the results obtained 
for clusters C^, Cjj, and Cj^j with the radius cores we 
have previously fixed, whereas the associated dotted lines 
correspond to the same clusters with smaller cores, the 
new core radius being 0.1 Mpc, 0.05 Mpc, and 0.025 AIpc 
in cases , C^^ , and Cj^^ , respectively (the initially cho- 
sen radius have been reduced by a factor 1/2). 



The continuous and dotted associated lines of Fig. |3| 
are very similar for any cluster, which implies that, for 
6 — 2.5', the lens deviation S has a very weak dependence 
on the core radius Tc. The direction 9 — 2.5' has been ar- 
bitrarily chosen in the 6 interval where the lens deforma- 
tions are relevant (see Fig. nevertheless, other angles 
covering this interval have been also studied. For exam- 
ple, we have investigated the directions corresponding to 
the maximum deviation in each of the cases reported in 
Fig. n and lens deviations also depends weakly on the 
core radius. The same occurs for any 9 value, except for 
very small angles, which are not relevant for the angular 
scales under consideration (see below). In Fig. 13 we can 
also see that the region where 'V±(j> is significantly con- 
tributing to the integral of Eq. (O has a comoving size 
of various Megaparsecs for any cluster of interest. It has 
been checked that the same conclusion is also valid for 
any relevant 9 values. 

A scale for the computation of angular correlations 
has been chosen. In order to do that, the CMB angular 
power spectrum (Ci coefficients) due to lensing by non- 
linear structures has been estimated using CMBFAST. 
Results lead to the conclusion that 6*2(0;) correlations 
from this lensing are important for a < 13.5' {£ > 800); 
hence we must work with some scale £ > 800. The scale 
£ = 1600 {a = 6.75') has been chosen to calculate cor- 
relations because this scale has the following properties: 
(a) for £ = 1600, the effect of linear inhomogeneities - 
which should be reduced as much as possible- is smaller 
than that of the case I — 800 (compare the bottom pan- 
els of Fig. EJ, (b) according to CMBFAST estimations, 
lensing from nonlinear structures gives a Cieoo coefficient 
which is smaller than the maximum Ce only by a factor 
~ 0.6, and (c) the value £ — 1600 is in the range of 
the multipoles to be observed by PLANCK (the most 
sensitive projected experiment for CMB anisotropy mea- 
surements). 

According to previous comments, a box of 128 AIpc 
have been assumed as the basic one and, then N-body 
simulations have been performed inside this box with a 
standard PM code llSll . which was used, tested, and de- 
scribed in detail in jl3|. The resolution has been succe- 
sively increased starting from a very poor one of 2 Mpc. 
Of course, such a low resolution leads to non-linear struc- 
tures whose amplitudes (sizes) are smaller (larger) than 
those of the true clusters. These simulations only traces 
a good spatial distribution of extended structures with 
total masses as those of the clusters (which would be- 
come clusters for appropriate resolutions). Resolutions 
of 1 Mpc and 0.5 Mpc have been also considered, but no 
higher resolutions have been used by the reasons given in 
Sees. V and VI. 



IV. DEPARTURES FROM GAUSSIANITY 

Suppose an observer which is located at a certain point 
with comoving spatial coordinates . The main ques- 
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tion is: which is the deviation, S{xp,n), observed from 
point Xp in the n direction? 

If 4>j:{t) is the FT of the potential cj), namely, 

<^(-' = J d^ke-^^^-'m ^ (16) 

and 5^: is the FT of the density contrast [p — Pmo)/ Pmo, 
Equation leads to the following relation in Fourier 
space: ip^ — BS^/k^, where B = —pmQ/a. 

Using elemental Fourier algebra and Eqs. ©, and 
(|7|). the following basic equations -giving the required 
deviation- are easily obtained: 

^^(^p'-) - / , (17) 

where 

Fj:{n) = [ " W{X)B{X)e-'^''"S^{X)dX , (18) 

and k± — k ~ {n ■ k)n. 

According to Eq. H17|l . each component of d{xp,n) is 
the FT -extended to all the space- of a component of 
vector k±Fj:{n)/k'^. The integral of the r.h.s. of Eq. lfTH|) 
must be estimated to get the function F-^{n) involved 
in the FT. The integration variable A can be seen as a 
generalized time coordinate. Equations 117|) and H18|) 
are general and they describe how the photons deviate 
-following null geodesies- in a realization of the full uni- 
verse; nevertheless, such a realization is not available in 
practice, and we are constrained to use these equations 
in a fictitious periodic universe, let us now discuss this 
fact and its consequences in detail. 

Structure evolution is simulated with a certain N-body 
code (see Sec. Ill), which involves a box and a certain 
resolution (a network in position space) . This code has a 
certain time step (or A step) and, consequently, it gives 
data at a set of time values; among these data, we are par- 
ticularly interested in the function 5-j:{\) which appears 
in Eq. (|18|l . This function is easily evaluated at each 
node, k, in Fourier space, and at each time step U, but 
not at arbitrary times. The question is: Can we calculate 
the integral ((TH|l using only the data corresponding to the 
time discretisation of the N-body code? Fortunately, the 
answer is positive. We have studied the function 5j^{X) for 
many k values and, thus, we have verified that the term 
J(fc,A) = W{\)B{\)5j:{\) -which appears in Eq. lfTS|l - 
can be very well approximated by a straight line between 
Xi and Ai+i, namely, between two successive time steps 
of the N-body simulation. The equation of this line is of 
the form: 

h{k,X)^ A,{k)\ + B,{k) (19) 

where quantities Ai and Bi are obtained from 5j:{Xi) and 
%(Ai+i); namely, from quantities given by the N-body 



in two successive steps. If Eq. (|19|) is substituted into 
Eq. (|18|) . the integral of the r.h.s. can be analytically 
calculated in the interval (ti, tj+i). The addition of these 
integrals for all the N-body intervals gives function F^{'n). 
After this integration, the FT H17|l can be performed and, 
thus, the deviation field in the n direction is simultane- 
ously calculated for many observers located at the nodes 
Xp of the spatial Fourier box. The importance of this fact 
is discussed below. 

Discretisation implies that the CMB photons do not 
move in a true realization of the universe, but in a pe- 
riodic universe, which is, at each time, an ensemble of 
boxes identical to that of the N-body simulation. This 
periodicity should lead to errors in our estimations. In 
order to understand the problem and its solution, let us 
establish a two-dimensional (2D) analogy. The squares 
of Fig. 21 (boxes in the 3D case) have a size of 128 Mpc. 
These squares contain small circles (clusters), and let us 
imagine an observer, S, at the center of the right square. 
If a photon propagates along the pointed line of Fig. 0] it 
crosses all the squares in the same way, namely, it passes 
at the same distances from the same clusters and under- 
goes the same systematic repeated deviations. Only evo- 
lution induces small deformations from square to square, 
but these deformations would be only relevant after a 
large enough number of squares have been crossed; hence, 
a certain wrong accumulative effect should appear. How- 
ever, if the photon moves following the dashed line of Fig. 
^ the situation is very different because the photon en- 
ters each square around a very different place and, then, 
it can be affected by distinct independent clusters. The 
key point is the selection of the angle $ and its associate 
distance L. Two facts are relevant to make this choice: 
(i) clusters significantly deviate the CMB photons when 
their impact parameter is smaller than a few Megaparsecs 
(see Fig. O, and (ii) correlations in the cluster distribu- 
tion (superclusters) extend near 10 Mpc; hence, if we 
choose L between 15 Mpc and 20 Mpc, photons move 
from square to square, but they are always influenced by 
different almost uncorrelated clusters and, consequently, 
the effect of periodicity must be minimum. All these con- 
siderations are easily generalized to the 3D case, in which 
there are two angles defining the line of sight. Hereafter, 
these angles are identified with the spherical coordinates 
and (j) (with respect to the Cartesian axis defined by the 
edges of the box) . These angles can be easily chosen to 
maximize the length of the path where the photons un- 
dergo the action of distinct uncorrelated clusters; namely, 
to minimize the effects of periodicity. In the 3D case, with 
boxes of 128 Mpc, it has been found that for 9 — 10.75° 
and (j) — 41.19°, the photons would cross 56 boxes before 
arriving to the initial region of the box. The CMB pho- 
tons moving along this direction (hereafter denoted Nq 
and called the preferred direction) do not feel periodicity 
during their travel from r ~ 7400 Mpc to r = 0, namely, 
from z — 3.62 io z — Q and, taking into account that 
the lens effect we are evaluating is produced at these low 
redshifts by galaxy clusters, it follows that, for direction 
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FIG. 4: Sketch of a 2D universe filled by squares which contains clusters (circles) and photon patches (lines). Dotted (dashed) 
line corresponds to a direction parallel (non parallel) to the square edges. 



iVo and close ones, periodicity must be rather irrelevant 
(see Sec. V for more discussion). Furthermore, the pho- 
ton enters in two successive boxes through points that 
-when placed in a unique box- are located at a distance 
L — 15.9 Mpc. For the chosen values of 9 and cj) and 
a box of 256 Mpc, we have L = 31.8 Mpc, and photons 
does not feel periodicity at all from decoupling to present 
time. 

On the angular scales under consideration, the clus- 
ter distribution producing lensing is statistically inde- 
pendent on the structure distribution causing the pri- 
mary anisotropy; in fact, for a values between 6.75' and 
13.5' (1600 > i > 800), primary anisotropy is due to 
temperature fluctuations previous to decoupling, which 
appeared as a result of the coupling between matter an 
radiation in the presence of energy density fluctuations; 
hence, the scales of temperature fluctuations and those 
of energy perturbations coincide. Inhomogeneities with a 
comoving size of 41 Mpc (20.5 Mpc) located at z = 1100 
subtend an angle of 13.5' (6.75'); hence, taking into ac- 
count that -in our ACDM model- the comoving scale of 
an Abell cluster with lO^^M© is 25.6 Mpc, we can con- 
clude that, for 800 < £ < 1600, the primary anisotropy is 
produced by comoving scales similar to those of clusters. 
These small structures located close to our last scatter- 
ing surface (z ~ 1100) and the clusters producing lensing 
(small z values) are very far, different, and independent 
structure distributions; therefore, we can consider that 
the distribution of clusters is producing deformations on 
independent primary CMB maps. In order to calculate 
the averages of Eq. I|15|) (correlations), various full real- 
izations of the cluster distribution should be considered 
to produce deviation fields and, then, each of these fields 
should be used to deform a large enough number of inde- 
pendent primary maps. In practice (computational lim- 
itations), only a few pairs of directions (close to A^o to 
minimize periodicity effects) are studied for each cluster 
distribution (N-body simulation), and these directions 
should be used to deform different primary maps. Let us 
now describe in detail the method implemented in this 



paper to calculate correlations. 

For each N-body simulation, nine directions Hi are con- 
sidered. These directions are chosen in such a way that 
they depict the vertices of four neighboring squares with 
size a on the last scattering surface. These squares form 
another bigger one with size 2a with vector Nq point- 
ing towards its center. It is important that, for each of 
the chosen directions, the deviations 5{xp,ni) are calcu- 
lated for the observers, Xp, at the same time (see 
above). Let us discuss the importance of this fact from 
the statistical point of view. Put a 4.25° x 4.25° map, 
Mp, of primary anisotropy in such a way that, from the 
box center (point x^), the direction A^o points towards 
the map center C . Consider now one of the observation 
direction n^. From the box center, the chosen direction 
points towards a certain point Ci, nevertheless, from an- 
other point of the box x^ , direction Hi does not point 
towards point Ci but towards another point Qi of the 
map Mp. Two parallel lines -with direction Hi- starting 
from points x^ and x^ of the box would intersect the 
map Mp in points Ci and Qi, respectively; hence, the Qi 
coordinates in the map Mp can be calculated as follows: 
project the vector Xq — x^ on the plane normal to Hi, 
translate this projection parallely from the hypersurface 
t = to the last scattering surface at z = 1100, calcu- 
late the angle subtended by the final vector (as it is seen 
by the observers) and, then, the angular coordinates of 
point Qi in the map Mp can be trivially obtained. Of 
course, the deviation, 6{x^,Hi), corresponding to the ob- 
server x^ in the direction Hi, must be applied to point 
Qi in the map Mp . Hence, for each direction Hi , we 
have not a unique point in Mp, but many points, one 
for each x^^ observer, and these points cover a certain 
region on map Mp . In order to estimate the size of the 
covered region, we can consider a fixed direction, for ex- 
ample H = (0, 0, 1), and a box of 128 Mpc (256 Mpc) per 
edge; thus, given two observers x^-^ and x^^ separated by 
a distance of 128 Mpc (256 Mpc), the angular distance 
between points Qi and Q2 in the map Mp appears to 
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be S = 0.5° (S = 1°), which means that, for a = 13.5', 
a squared region with a size of 43.5' (73.5') is covered, 
and for a — 6.75', the size of the corresponding region is 
36.75' (66.75'); hence, we see that, for boxes of 128 Mpc 
(256 Mpc) and a = 6.75', the size of this region is larger 
than a by a factor close to five (ten) and, consequently, 
the averages in the observers ~a novel aspect of our 
ray-tracing procedure- play a crucial statistical role. The 
greater the boxes and the smaller the angular scales a, 
the larger the above factor and the greater the statistical 
significance of the averages. 

On account of the above considerations, correlations 
are calculated as follows: (i) twelve different pairs of di- 
rections, sixteen triads, and four tetrads are defined us- 
ing the above nine directions, (ii) given an observer x^ , 
a map (primary anisotropy), and a N-body simu- 
lation (for S computation), the deviations S{xQ,ni) and 
Eq. 1^ (or Eq. 10)) are used to calculate the lens ef- 
fect A[ni + (5(x^,nj)] - A[ni] (or 6{x^,n{) ■ [8/^^/ dn]fi.) 
for each of the Hi directions and, then, the resulting data 
are used to compute the averages H15II for the above pairs 
(to = 2), triads (to = 3), and tetrads (m — 4), (iii) cal- 
culations in (ii) are repeated for each observer Xq , and 
results are averaged again, (iv) all the process is repeated 
for a certain number of Mp maps and a new average is 
performed and, (v) finally, the above calculations are re- 
peated for a certain number of N-body simulations, and 
the final averages are taken as our estimations of the cor- 
relations 115|l . The question is: How many primary maps 
and how many cluster realizations would be necessary to 
calculate the required correlations? 

In order to answer this question, the following method 
has been implemented: In a first step, the process (i)- 
(iv) has been used to analyze pure primary maps in the 
absence of lens deviations. The primary maps are not de- 
formed and, consequently, the resulting averages should 
approach the true two, three, and four direction correla- 
tions of these maps, which are hereafter called primary 
correlations. These correlations can be estimated by a 
direct analysis of the simulated primary maps (see next 
section) and the results of such an analysis can be com- 
pared with the correlations estimated by the process (i)- 
(iv) in the absence of deviations. Since no deviations 
have been considered, this study measures the capability 
of our method (based on nine direction, observers, 
and various primary maps) to create pairs, triads and 
tetrads covering a significant part of the primary maps; 
namely, allowing a good statistical analysis of maps. In 
the second and final step, the number of primary maps 
suggested by our previous study (in the absence of devia- 
tions) is taken and, then, deviations from more and more 
N-body simulations are considered in order to compute 
the required correlations as it has been described above. 
When the averages reach almost stable values, the pro- 
cess is stopped and no more N-body realizations are per- 
formed. Results are presented in next section. 



V. RESULTS 

First of all, the methods applied -in this paper- to cal- 
culate correlations of the lens effect are tested. In order 
to do that, these methods are used to analyze well known 
maps of primary anisotropy. The angular scale for cor- 
relation estimates is a = 6.75'. A set of two hundred 
maps of primary anisotropy has been used. Results are 
presented in Fig. |S1 where variable Np (horizontal axis) 
is the number of primary maps used to get the correla- 
tion value appearing in the vertical axis. Left, central, 
and right panels correspond to the two (m = 2), three 
(to — 3) and four (to = 4) direction correlations l|15|l . 
respectively. In all the panels, the solid line shows the 
correlations obtained by means of an exhaustive coverage 
of the Mp maps using many pairs, triads and tetrads of 
directions. This is our best estimate of the required cor- 
relations; results from this method -which cannot be ap- 
plied to analyze the lens component- approach the theo- 
retical values of the correlations we have used to simulate 
the Mp maps. Other methods must give similar correla- 
tions to be acceptable. As Np increases, the solid lines 
of the top and bottom panels seem to tend to a certain 
value, whereas the continuous line of the middle panel 
seems to be compatible with a vanishing C3 correlation. 
The long dashed lines display the correlations obtained 
using nine directions and only one observer (only nine di- 
rections in each map). These lines slowly approaches the 
solid one as Np increases and the level of approximation 
seems to be good for Np = 200. The dotted (dashed) 
lines give the correlations obtained with nine directions 
and 128'^ (256'^) observers. As it follows from Fig. [SJ 
dashed and dotted curves approach the solid ones faster 
than the long dashed lines and, consequently, the use of 
many observers located in the nodes of a network (see 
Sec. IV) is statistically significant; the greater the box 
size, the faster the statistical approximation to the solid 
lines. The best situation corresponds to 256'^ observers 
in a 256 AIpc box. Two hundred primary maps suffice 
for a good enough estimate of C2(6.75) and 6*4(6.75); 
particularly, when N^ observers are considered. 

After computing the correlations C2(6.75), C3(6.75), 
and 6*4(6.75) for lens deformations, it would be worth- 
while a comparison between the correlation level of 
these deformations and that of primary maps. The ra- 
tios r2 = C2(6.75)/C72(0), = C3(6.75)/C|''^(0), and 
r4 = C4(6.75)/C|(0) are apropriate to do this compar- 
ison; these quantities will be calculated for lens defor- 
mations and primary maps and, then, results will be 
compared. We begin with the case of primary maps. 
Using two hundred of these maps and exhaustive cov- 
erage, we have found: Ci = {Ap{n)) = —2.12 x 10^^, 
62(0) = {Ap{n)Ap{n)) = 7.29 x 10-^ and n = 

Ci/Cl^^{0) = 7.85 X IQ-'^. We know that there are 
many positive and negative temperature contrasts -in 

1 /2 

the maps- with absolute values close to (0) and, 
consequently, only an strong cancellation of these con- 
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FIG. 5: Correlations C2, C3, and C4 extracted from simulated maps of the primary CMB anisotropy (for a = 6.75'), versus 
the number Np of simulated maps used in the computations. Left, central, and right panels show the correlations C2, C3, and 
C4, respectively. Different curves correspond to distinct methods for correlation estimates which are described in the text. 



trasts can explain the resulting value of the mean con- 

1/2 

trast Ci(0), which is much smaller than (0) (small 
ratio ri ^ 7.9 x 10"^). Strong cancellation leading to 
a small ri value indicates that, as the number of maps 
increases, quantity {Ap{n)) approaches either a small 
value or zero. Similarly, in order to calculate the correla- 
tion C2(6.75) = (A^(?Tri)Ap(7i2)), positive and negative 
numbers with absolute values close to 6*2(0) are aver- 
aged. The resulting ratio 7-2 ~ 0.85 indicates that these 
numbers have not cancelled among them and a nonva- 
nishing 6*2(6.75) correlation exists. For m = 3 (m=4), 
the relevant ratio is ^ 8.8 x 10~^ (r4 ~ 2) and, this 
small (large) ratio 7-3 {r^) strongly suggests a vanishing 
C3 (significant C4) correlation. 

Since the simulated primary maps are Gaussian, the 
relations 6*3(6.75) = and 6*4(6.75) = 36*1(6.75) must 
be satisfied. The correlation 6*3 obtained by exhaustive 
coverage (solid hues of Fig. [SJl is 6*3 ~ —1.76 x 10~^^, 
and it has been proved to be small by computing r^; fur- 
thermore, we have used the same method to get the cor- 
relations 6*2 ~ 6.24 X 10-8 and 6*4 ~ 1.09 x 10" i"*, which 
lead to the ratio 36*|/C4 ~ 1.07; hence, we can say that 
the correlations extracted from the maps are compatible 
with Gaussianity, as it should be for the primary maps 
we are analyzing. If we use other methods to analyze 
the Mp maps; for example, the methods used to built 
up the dotted, dashed, and long dashed curves of Fig. O 
(see above), the extracted correlations appear to be sim- 
ilar to those reported above and so they are compatible 
with Gaussianity; this means that all the methods used 
to analyze maps work very well for primary maps, and 
the same should occur in our applications to the analysis 
of lens deformations. 

Correlations of lens deformations have been calculated 
in various cases using distinct PM simulations; in case 1, 



the resolution is 0.5 Mpc (the best one), the box size is 
128 Mpc, directions are close to the preferred one, and 
the angular scale is a = 6.75'. Correlations obtained in 
this case are presented in the first row of table P and 
also in the continuous lines of the left panels of Fig. 
All the data of table are mean correlations with 2a 
errors calculated from the correlations of forty simula- 
tions. These errors are statistical ones assigned in case 
1, but there are other sources of uncertainty as period- 
icity and simulation features, which are not included in 
the errors reported in table The lines mentioned above 
display the variation of the estimated correlations as the 
number N of simulations (used to perform averages) in- 
creases; for N > 20, correlations 6*2 and 6*4 seem to 
coverge towards a certain value, whereas the correlation 
6*3 slowly decreases in the iV-interval (20,40) (see con- 
tinuous lines in the middle panels of Fig. El which are 
two different representations of the same function). We 
assign to 6*2, 6*3, and 6*4 the values and errors given in 
the first row of table U the tabulated values of 6*2 and 
6*4 are estimations of these correlations, whereas that of 
6*3 can be better considered as an upper limit due to the 
slow decreasing we have mentioned above. The values 
assigned to 6*2, 6*3, and 6*4 have been used to calculate 
the ratios ~ 0.081, rs ~ 1.17 x 10"^, and r4 ~ 0.014, 
which are smaller than the corresponding ratios of the 
primary maps (see above). The small value indicates 
that lensing (by nonlinear structures) does not introduce 
any appreciable 6*3 correlation, whereas the small values 
of 6*2 and 6*4 suggest small levels of these correlations. 
Finally, the ratio 36*|/6*4 ~ 1.43 confirms a deviation 
with respect to Gaussianity. 

Since the method proposed in previous sections is 
based on N-body simulations, results could depend on 
both resolution and box size. These dependences are 
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FIG. 6: Correlations C2, C3, and C4 extracted from simulated maps of the lens effect (for a = 6.75'), versus the number N of 
simulated maps used in the calculations. Top, middle, and bottom panels show the correlations C2, C3, and C4, respectively. 

solid, dashed, and dotted lines of the left panels corresponds to resolutions of 0.5 Mpc, 1 Mpc, and 2 Mpc, respectively. Dashed 
(solid) lines of the right panels show correlations for direction close to the preferred one (a direction parallel to an edge of the 
simulation box). 
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TABLE I: Correlations in simulated lens deformations of the CMB sky. 



Case 


Box Size 
in Mpc 


Resolution 
in Mpc 


Direction 


C2 X 10" 


Ci X 10'« 


Ci X 10^' 


1 


128 


0.5 


optimal 


3.59 ± 0.12 


-1.10 ± 0.14 


2.71 ± 0.18 


2 


128 


1.0 


optimal 


3.57 ± 0.06 


-1.40 ± 0.18 


2.70 ± 0.10 


3 


128 


2.0 


optimal 


3.27 ± 0.06 


-2.63 ± 0.24 


2.28 ± 0.08 


4 


128 


1.0 


parallel to an edge 


2.50 ± 0.02 


-1.43 ± 0.16 


1.05 ± 0.03 


5 


256 


1.0 


optimal 


3.08 ± 0.06 


-0.76 ± 0.06 


2.07 ± 0.08 



TABLE II: Comparing lens correlations for various pairs of 
estimations. 



Compared Cases 


IOOAC2/C2 


IOOAC3/C3 


IOOAC4/C4 


1-2 


0.56 


21. 


0.37 


2-3 


8.91 


58. 


16. 


2-4 


30 


2.09 


61 


2-5 


13.73 


46 


23 



now studied. In order to analyze the importance of res- 
olution, cases 2 and 3 have been considered, they are 
identical to the case 1 described above except for res- 
olution, which is 0.5 Mpc in case 1, 1. Mpc in case 2, 
and 2. Mpc in case 3. Results of case 2 (3) are dis- 
played in the second (third) row of table and also in 
the dashed (dotted) lines of the left panels of Fig. El 
where we see that the correlations C2 , C3 , and C4 of case 
3 (2. Mpc resolution) are smaller than those of case 2 
(1. Mpc) and case 1 (0.5 Mpc), which are very similar 
between them. In order to quantify these comparisons 
between results corresponding to pairs of cases, the rela- 
tive variations ACm/C*™ for m = 2,3,4 are obtained for 
the pairs 1-2 and 2-3, and the comparation percentages 
100AC„i/C,„ for these pairs are given in the first and 
second rows of table |nl where we see that the compara- 
tion percentages of the pair 1-2 are much smaller than 
those of 2-3. These considerations strongly suggests that 
a 2. Mpc resolution is poor, whereas a 1. Mpc resolution 
is good because a better one (0.5 Mpc) does not lead 
to significantly better results; particularly, for the cor- 
relations C2 (IOOAC2/C2 = 0.56 for cases 1-2) and C4 
(IOOAC4/C4 — 0.37 for cases 1-2). Taking into account 
that cluster description can be improved using higher res- 
olutions, these results have two possible interpretations: 
(1) the calculation of the lens correlations does not re- 
quire a more detailed description of the clusters (for the 
chosen angular scale a = 6.75') and, (2) results from 
resolutions of 1 Mpc and 0.5 Mpc are similar because 
of a lack of resolution in both cases, but higher resolu- 
tions would lead to -different- better results. We have 
not found any theoretical reason supporting one of these 
alternatives. In this situation, the most effective proce- 
dure is the use of new simulations with higher resolutions, 
which should support one of the above possibilities (see 
Sec. VI for more discussion). 

After concluding that a resolution of 1. Mpc seems to 



be good (at least it is equivalent to a better resolution 
of 0.5 Mpc), and before analyzing the importance of the 
box size, a study of the relevance of periodicity is worth- 
while. In order to perform such a study, results of case 
2 are compared with those of case 4. In both cases, the 
same forty simulations are used, but the propagation di- 
rections of the CMB photons are different. In cases 2 and 
4, directions form small angles with the preferred direc- 
tion and with an edge of the simulation box, respectively. 
Results of case 4 are displayed in the forth row of tabled 
and also in the solid line of the right panels of Fig. |H| 
(where dashed lines correspond to case 2). Comparing 
the solid and dashed lines of these panels, one easily con- 
cludes that periodicity strongly (weakly) affects the cor- 
relations C2 and C4 (C3). Accordingly, the comparation 
percentages of the correlations C2 and C4 for the pair of 
cases 2-4 (third row of table |nj are much greater than 
that of the C3 correlation, which appears to be as small 
as 2.09%. This dependence of C2 and C4 on the direction 
constraints us to work with directions close to the pre- 
ferred one, which strongly minimizes periodicity effects 
(see Sec. IV). 

Finally, results of case 2 are compared with those of 
case 5 to analyze the importance of the box size. Case 5 
is identical to case 2 except for this size, which takes on 
the values 256 Mpc and 128 Mpc in cases 5 and 2, re- 
spectively. Possible differences between cases 2 and 5 can 
be explained in two different ways: (i) boxes of 256 Mpc 
contain linear inhomogeneities more extended than those 
involved in boxes of 128 Mpc and, according to our dis- 
cussion of Sec. Ill, linear inhomogeneities with different 
sizes produce distinct lensing and, (ii) there is a resid- 
ual effect of periodicity, which would depend on box size 
for the fixed preferred direction used in our calculations 
(initially associated to 128 A-Ipc boxes, see Sec. IV). Re- 
sults of case 5 are shown in the fifth row of table and 
also in the solid lines of Fig. [7| (where dashed lines corre- 
spond to case 2). Comparing the solid and dashed lines 
of Fig. and taking into account that the comparation 
percentages of the pair 2-5 (forth row of table are 
13.73% for C2, 46% for C3, and 23% for C4, we conclude 
that box size is important. Is that due to residual peri- 
odicity? The answer to this question is negative because 
the comparation percentage in C3 due to periodicity has 
been proved to be very small (2.09%), whereas the cor- 
responding percentage has appeared to be 46% when the 
box size has been varied. Hence, the box size is important 
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FIG. 7: Same as Fig. |^for box sizes of 128 Mpc (dashed lines) and 256 Mpc (solid lines). 



as a result of the presence of large scale inhomogeneities 
and, consequently, our best estimate of the correlations 
induced by nonlinear lensing corresponds to the smallest 
box size (128 Mpc) and the best resolution (0.5 Mpc) we 
have used (case 1). The comparation percentages of the 
pair 2-5 are not to be confused with the error percentages 
in the correlations of case 2, which would be expected to 
be smaller, perhaps about a few per cent. The same has 
occurred for resolution, where the comparation percent- 
ages of the pair 2-3 (second row of table In]) are much 
larger than the error percentages of case 2 (comparable 
to the quantities in the first row of table lll)l . 



VI. DISCUSSION 

A ray-tracing method to calculate correlations in lens 
deformations of CMB maps has been designed, tested, 
and applied. It is based on the use of PM N-body sim- 
ulations to evolve nonlinear cosmological structures. Al- 
though correlations are calculated from the deformations 
associated to a reduced number of directions (no com- 
plete maps of deformations are created) , numerical calcu- 
lations are time consuming because a large enough num- 
ber of N-body simulations must be performed to calcu- 
late good statistical averages (correlations) . The method 
is similar to that proposed by Aliaga et al. 3 to study 
the Rees-Sciama effect, and it can be easily implemented 
starting from an N-body code. 

Firts of all, the main structures producing lens defor- 
mations in the CMB sky have been identified. A qual- 
itative estimation of the lens deformations produced by 
a representative set of linear and nonlinear cosmological 
structures has been used to achieve this identification. 
Various reasons have motivated the choice of the angular 
scale a = 6.75' to compute correlations (sec Sec. III). 
For this scale, it has been proved that lens deformations 
are mainly produced by big and standard clusters, al- 
though there is a moderate contribution from linear inho- 
mogeneities with diameters greater than 100 Mpc. This 
information has been crucial to choose the N-body simu- 
lations in such a way that the effect of the linear inhomo- 



geneities contained in the simulation box is much smaller 
than that of the nonlinear structures. Using these simu- 
lations, the deviation field 5 produced by nonlinear struc- 
tures can be numerically estimated and, independently, 
the field 5 due to large scale linear inhomogeneities can 
be analytically calculated. Qualitative arguments sug- 
gest that box sizes between 128 Mpc and 256 Mpc could 
be appropriate to achieve the above requirements (Sec. 
Ill) and, then, experiments with simulations (see Sec. V) 
have been used to reach the following conclusions: (i) a 
size of 128 Mpc is preferable against larger sizes including 
extended linear structures (which would produce a signif- 
icant lensing) and, (ii) resolutions of 0.5 Mpc and 1. Mpc 
lead to similar correlations; hence, results converge as res- 
olution increases. Is this convergence real or apparent? 
As it has been pointed out in Sec. V, the resolution of our 
simulations must be increased to answer this question. 
How much should we increase the resolution? We should 
consider -at least- a resolution of 0.25 Mpc and, per- 
haps another of 0.125 Mpc. These resolutions lead to se- 
vere computational problems because many simulations 
are necessary for statistical purposes. Higher-resolution 
codes -for example tree or P^M codes- could be more ap- 
propriate than our PM one. Fortunately, our simulations 
lead to boxes containing nonlinear structures which allow 
us to apply and test our ray-tracing procedure and, fur- 
thermore, our best resolution (0.5 Mpc) should lead to a 
rather acceptable description of clusters, because the cen- 
tral region, where the density is more peaked (distances 
^ 1 Mpc from the center), is covered by about 2^ pairs 
of cells. We can state that our main goal (description 
of a new ray-tracing procedure and estimations of non- 
gaussianity from nonlinear structures) has been reached, 
whereas results from higher resolutions are out of the 
scope of this work, and they will be developed in due 
time and presented elsewhere. 

In our model, photons move in a periodic universe 
which is the repetition of the simulation box. We have 
verified that periodicity is important, but photons do not 
feel periodicity when they move along directions close 
to a preferred one (see Sec. IV); this fact constraints 
us to calculate deformations for this type of directions. 
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Fortunately, these directions suffice to find the required 
correlations. It is worthwhile to emphasize some novel 
aspects of our calculations: the basic equations are ifTTjl 
and (|18|) : the r.h.s of Eq. p7|) has the form of a Fourier 
transform, and it gives -at the same time- the deviation 
field S{xp, n) for all the observers located in the nodes ( 
Xp ) of the simulation grid. This fact has been probed 
to be relevant when averages are performed. The r.h.s. 
of Eq. (|18|l is a time integral which can be analytically 
calculated after using the interpolation formula (|19|l , this 
procedure simplifies the calculations. The use of direc- 
tions close to the preferred one is the last novel aspect 
deserving attention. For these directions and fixed boxes, 
photons enter the n-th box across a certain region and 
the (n-l-l)-th box through a different independent zone; 
hence, roto-traslations are not necessary to avoid period- 
icity effects; furthermore, it is worthwhile to emphasize 
that -without roto-traslations- function V^i/) is contin- 
uous at the points where photons cross boundary boxes 
(periodicity). However, if the boxes are moved (roto- 
tralations), discontinuities at these points are unavoid- 
able and, consequently, the function to be integrated in 
Eq. (5) has artificial finite discontinuities which influence 
deviation calculations. Perhaps these discontinuities do 
not produce any relevant total effect when many boxes 
are crossed, but potential problems with discontinuities 
are clearly surmounted by our method in the absence of 
roto-traslations, which seem not to be either necessary 
or appropriate for us. 

Since periodicity effects are prevented along any direc- 
tion forming an angle of a few degrees with the preferred 
one, our method can be used: (1) to simulate squared 
maps with sizes of a few degrees, for example, 4° x 4° 
maps with the preferred direction pointing towards their 
centers, and (2) to compute angular correlations C„ for 
n > 4. In previous sections, correlations C3 and C4 have 
been computed using sets of three and four directions 
close to the preferred one and, evidently, using the same 
methods and other sets of these directions, the correla- 
tion Cg and higher ones can be computed. On account of 
these considerations, future applications of the proposed 
method are certainly promising. 

Hereafter, we only make reference to the lens effect 
(or lens deformations) produced by nonlinear structures. 
As it was predicted by Bernardeau 0, the lens effect 
does not create any signiflcant C3 correlation (ra ~ 
1.17x10"'^). Correlations C2 and C4 appear to be weaker 
than those of the primary maps (low correlation level, see 
Sec. V). The resulting ratio 3C|/C4 ~ 1.43 points out a 
deviation from Gaussianity. In the observational maps, 



the lens deformations are superimposed to both the dom- 
inant primary anisotropies and other contributions de- 
scribed in the introduction (foregrounds and nonlinear ef- 
fects). Let us analyze the superposition of primary maps 
and lens deformations taking into account that: (1) the 
correlation C2 (^4) of primary maps has appeared to be 
greater than that of lens deformations by a factor close 
to 1700 (4 X 10^), (2) the correlation C3 vanishes in pri- 
mary maps and it has been found to be negligible for 
deformations, and (3) the two components under consid- 
eration are statistically independent. Simple calculations 
give the correlations of the superposition of components, 
which are almost identical to those of the primary maps. 
Deviations with respect to Gaussianity are almost neg- 
ligible in this superposition. In the observational maps, 
there are other deviations due to the presence of other 
components, which would contribute to hide the small 
deviations associated to lens deformations. 

In future, our calculations could be repeated using 
other N-body codes and resolutions, the lens effect due to 
linear structures contained in boxes with different sizes 
could be studied in more detail; the present study could 
be extended to other components of the observational 
CMB maps, in particular to the Rees Sciama compo- 
nent |3|; maps of lens deformations and other compo- 
nents could be simulated -with our method- and these 
maps (which would have the true statistical properties at 
any order of correlation) could be analyzed using various 
estimators of nonGaussianity proposed in the literature 
ph. .21] ; our simulations could be also used to study the 
correlation between galaxy shear and CMB temperature 
distortions 0| and, finally, using other N-body codes 
and greater resolutions, it would be worthwhile the study 
of very small (large) angular scales (£ values), for which, 
the Ce coefficients of the primary anisotropy are very 
small. Some of these studies and other possible ones 
could be useful to improve on the method proposed in 
this paper; nevertheless, a good version has been already 
implemented and numerical codes are operative. 
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